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Quantum magnetism lies at the heart of 
many intriguing phenomena in condensed mat- 
ter physics. Its manifestations range from an- 
tiferromagnets to spin-hquids, and it is be- 
heved to play a central role in high-temperature 
superconductivity^ ^. Remarkably, even simple 
models of the underlying many-body physics are 
often intractable with current theoretical and 
computational methods. The controlled setting 
of ultracold fermionic atoms in optical lattices is 
therefore regarded as a promising route to pro- 
vide new insights^ ^. Yet, the low temperature 
scale required for entering the regime of quantum 
magnetism has hindered progress for optical lat- 
tice based systems. So far, superexchange oscil- 
lations on isolated double wells, one-dimensional 
decoupled Ising spin chains and classical mag- 
netism on a triangular lattice were studied with 
bosonic quantum gases^ ^. Here we report on the 
observation of quantum magnetism of a Fermi 
gas in an optical lattice. The key to obtaining 
and detecting the short-range magnetic order is 
a tunable geometry optical lattice. When load- 
ing a low-temperature two-component gas with 
repulsive interactions into either a dimerized or 
anisotropic simple cubic lattice, we find magnetic 
correlations on neighbouring sites. The correla- 
tions manifest as an excess number of singlets 
as compared to triplets consisting of two atoms 
with opposite spins. For the anisotropic lattice, 
we determine the transverse spin correlator from 
the singlet-triplet imbalance and observe antifer- 
romagnetic correlations along one spatial axis. 

Ultracold Fermi gases with repulsive interactions in op- 
tical lattices are an almost ideal realization of the Fermi- 
Hubbard model^^, having the potential to emulate its 
elusive low-temperature phase diagram. While density 
ordering has been explor ed exp erimentally in the metal- 
Mott insulator transition^^EU, spin ordering has so far 
not been observed owing to the low temperatures re- 
quired. The physics of the low-energy Hubbard sector 
is determined by the interplay of kinetic energy, inter- 
particle interactions and spin ordering, which is essen- 
tial to antiferromagnetic order, resonating valence-bond 
states, pseudo-gap physics and spin liquid^^^t^^. 

In this work, we use a tunable geometry optical 
lattic^^^^ to both access the low energy sector of the 
Hubbard model and probe magnetic correlations. A cen- 
tral requirement for quantum magnetism are tempera- 
tures below the exchange energy associated with spin or- 



dering. The atoms are either prepared in a dimerized or 
an anisotropic simple cubic lattice, as shown in Fig. [l] 
In both geometries, a subset of links of the underlying 
simple cubic lattice is set to a larger exchange energy as 
compared to the other links. This determines the charac- 
ter of magnetic ordei^^ and, given a fixed total entropy, 
allows for temperatures between the two exchange energy 
scales. In the dimerized lattice the resulting correlations 
on the strong links correspond to an excess population of 
the low energy singlet as compared to the triplet state. 
In the anisotropic lattice the low temperatures lead to 
antiferromagnetic spin correlations, the transverse com- 
ponent of which is also detected via a singlet-triplet im- 
balance. 

The experiment is performed with a harmonically con- 
fined, balanced two-component mixture of a quantum 
degenerate Fermi gas of ^^K. The atoms are prepared 
in two magnetic sublevels tuf = —9/2,-7/2 of the 
F = 9/2 hyperfine manifold, denoted by t and ^, at 
temperatures below 10% of the Fermi temperature. We 
load 50, 000 — 100, 000 repulsively interacting atoms at 
an 5- wave scattering length of 106(1) ao into the three- 
dimensional optical lattice, where ao denotes the Bohr ra- 
dius. The lattice is created by a combination of interfer- 
ing laser beams operating at 1064 nm with lattice depths 
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FIG. 1. Magnetic spin correlations. Schematic view of 
the nearest-neighbour spin correlations observed in the exper- 
iment. A two-component mixture of fermionic atoms (red and 
blue) is prepared close to half-filling in a cubic lattice with two 
different tunnel coupling configurations, a, Dimerized lattice 
with the strong dimer links td and weaker links t. Low tem- 
peratures lead to an excess number of singlets over triplets, 
b, Anisotropic lattice with strong and weak tunnelling ts and 
t along different spatial axes. Antiferromagnetic spin correla- 
tions in the transverse direction are formed along the strong 
link direction. In both figures exemplary thermal excitations 
in the form of spin excitations or holes are shown. 
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FIG. 2. Detection scheme. Summary of the technique used for measuring the atomic singlet and triplet fractions ps and pto- 
a, Schematic view of the different detection steps for the exemplary case of two singlet states in a dimerized lattice. Depending 
on the oscillation time, the absorption images on the right show either a large double occupancy corresponding to singlets (top 
row), or an increased higher band fraction indicating triplet states (bottom row), b. The two possible merging configurations 
in a dimerized lattice. Singlets and triplets are detected on a set of adjacent sites arranged on a chequerboard pattern in the 
plane, c. Exemplary singlet-triplet oscillation in a strongly dimerized lattice with U/t = 11.0(8) and td/t = 22(2). We observe 
an oscillation in the double occupancy after merging, ?7d, and in the higher band fraction, ?7hb, when merging along the strong 
links (blue data), whereas no oscillations are visible for the weak links (ochre data). The phase of the oscillation is shifted 
owing to the double occupancy removal procedure (Methods). The red dashed lines denote the extracted singlet and triplet 
fraction ps and ptQ. Error bars show the standard deviation of at least five measurements. 



Vx, Vy and Vz (Methods). We independently control 
the tunnelling strengths along all three spatial axes. In 
addition we can introduce a chequerboard dimerization 
in the xy plane by strengthening every second tunnelling 
link along the x axis, see Fig. [T^. The chequerboard 
pattern replicates along the z axis. Our system is well 
described by a three-dimensional single-band Hubbard 
model with repulsive on-site interaction energy /7, unless 
explicitly stated. The tunnelling along the weak links in 
both lattice geometries is set to t/h = 67(3) Hz for all 
measurements, where h denotes Planck's constant. 

As shown in Fig. [2^ and b, the fraction of atoms form- 
ing singlets and triplets on neighbouring lattice sites {p^ 
and pto) is detected by transforming the lattice to a 
chequerboard geometry, similar to a previously devel- 
oped technique^. In the dimerized lattice our detection 
scheme locally projects onto the two-site eigenstates of 
isolated dimers, whereas for the anisotropic lattice the 



singlet state is directly projected onto (It,!) — \iA))/V^ 
(Methods). In the first detection step, the atomic mo- 
tion in the initial lattice is frozen by rapid conversion 
to a simple cubic structure with negligible tunnelling. 
Next, all atoms on doubly occupied sites are removed. 
We then apply a magnetic field gradient, which creates 
a differential bias energy A for atoms of opposite spins 
on adjacent sites and causes coherent oscillations be- 
tween the singlet |s) = (|t,|) - ||,t))/V2 and the triplet 
|to) = (It, I) + state at a frequency u = A/h. 

If the initial amount of singlets and triplets is equal, no 
overall oscillation will be visible, as |s) and |to) oscillate 
in antiphase. 

After a certain oscillation time, we remove the gradient 
and merge two adjacent sites. Owing to the symmetry of 
the two-particle wavefunction, the singlet state on neigh- 
bouring sites evolves to a doubly occupied site with both 
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atoms in the lowest band, while the triplet state trans- 
forms into a state with one atom in the lowest and one 
atom in the first excited band. The fraction of atoms on 
doubly occupied sites t^d in the merged lattice is detected 
by a radiofrequency transfer to the previously unpopu- 
lated rriF = —5/2 spin stat^^. The fraction of atoms 
in the higher band 77HB is obtained from a band map- 
ping technique^. For the final readout we take absorption 
images after Stern-Gerlach separation of the spin states 
during ballistic expansion. For an imbalance between 
the initial singlet and triplet populations, 77D and 77HB 
will show oscillations with opposite amplitudes. As the 
double occupancy contains only contributions from two 
particles with opposite spins, we can infer the fraction 
of atoms forming singlets and triplets from the maxima 
and minima of a sinusoidal fit to 770- The higher band 
fraction has an additional offset caused by dimers con- 
taining two atoms with the same spin or one atom with 
an antisymmetric spatial wavefunction. 

When loading atoms into a strongly dimerized lattice 
and merging along the strong links, we observe oscilla- 
tions in T^D and t^hb, see Fig. This reveals an excess 
number of singlets, corresponding to magnetic order on 
neighbouring sites. We quantify this order by the nor- 
malized imbalance 

j-^ Ps-Pto ^ (1) 
Ps + Pto 

The order in the strongly dimerized lattice originates 
from temperatures below the intra-dimer exchange en- 
ergy Jd = -U/2 + y/lQtl + [72/2, which denotes the 
singlet-triplet splitting on a single dimer. While such 
temperatures require very low entropies for isotropic 
latticeJ^, in our system the access to the regime of 
magnetic ordering is facilitated by the presence of the 
weaker inter-dimer exchange energy J <C Jd- This leads 
to an entropy redistribution from states on the strong 
to the weak links and gives access to the temperature 
regime J < k^T < for experimentally attainable en- 
tropies (here denotes the Boltzmann constant). As 
expected for strong dimerization, we find no visible oscil- 
lations when merging along the weak links, which indi- 
cates the absence of magnetic correlations on these links, 
see Fig.[2]3. The observed constant values of t^d = 0.07(1) 
and 77HB = 0.63(3) are consistent with a state where 
nearly all singlets are located on neighbouring dimer 
links, with vanishing correlations between them. 

To analyse the effect of temperature on the magnetic 
correlations, we measure the dependence of the singlet- 
triplet imbalance on entropy, see Fig. [3^. The im- 
balance X and the absolute singlet fraction reduce 
for larger entropies, as triplet states become thermally 
populated. The singlet fraction is additionally dimin- 
ished by a shrinking half-filled region in the trapped 
system^^. We find good agreement with a second or- 
der high-temperature series expansion of coupled dimers 
when including an entropy increase of As = 0.4 /cb with 
respect to the initial entropy in the harmonic trap, s^^. 
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FIG. 3. Dimerized simple cubic lattice, a, Singlet-triplet 
imbalance on the strong dimer links vs. initial entropy be- 
fore loading into the lattice Sin and temperature k^T / Jd in a 
dimerized lattice with U/t = 11.0(8) and U/t = 22(2). The 
imbalance and the absolute singlet fraction (inset) decrease 
with increasing entropy. Solid lines are the prediction of a 
high-temperature series expansion taking into account differ- 
ent amounts of added entropy As during the lattice loading, 
b, Imbalance versus dimerization td/t and Jd/J, showing an 
increase for strongly dimerized simple cubic lattices. The solid 
line is the theory prediction for an entropy per particle of 
1.8 /cb in the lattice, which includes the heating during load- 
ing. Vertical error bars denote the fit error from singlet-triplet 
oscillations consisting of 63 measurements, the errors in td/t 
stem from lattice calibration uncertainties and the errors in 
Sin are the standard deviation of five measurements. Individ- 
ual curves of ps and pto for all measurements can be found in 
the Supplementary Information. 



This heating is associated to the lattice loading^^ and 
is larger for the lowest entropies, consistent with previ- 
ous results^^. From the measured imbalances we infer 
temperatures below 0.47^. 

For reduced dimerizations the coupling between dimers 
leads to increased inter-dimer correlations. The excita- 
tion energy of triplets is then lowered as they delocalize 
over the lattice, thus changing the nature of magnetic 
ordering. In Fig. we use the tunable lattice to in- 
vestigate the dependence of the imbalance X on the tun- 
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FIG. 4. Nearest-neighbour antiferromagnetic order. 

a, Transverse spin correlator versus tunnelling ratio tg/t and 
lattice depth in a three-dimensional anisotropic simple cu- 
bic lattice with Vy,z = 11.0(3) ^r. Positive values correspond 
to antiferromagnetic ordering. The inset shows the normal- 
ized spin correlator <S, denoting the fraction of antiferromag- 
netic ordering at the relevant density. Here U/t decreases 
from 16(1) to 10.5(8). Solid lines are the prediction of a high- 
temperature series expansion for an entropy per particle of 
1.8 /cb, as used in Fig.jsJ), and are shown up to ts/k^T — 1/2. 

b, Transverse spin correlator versus entropy before loading 
into the lattice dX U/t — 10.5(8) and ts/t = 7.3(6), together 
with a schematic view of the spin ordering. Error bars as in 
Fig.H 



nelling ratio td/t. As the dimerization is progressively 
removed the imbalance decreases in good agreement with 
theory and eventually falls below our experimental reso- 
lution. This decrease can be attributed to the inter-dimer 
exchange energy becoming smaller than the temper- 
atur T. For vanishing temperatures the system is ex- 
pected to undergo a quantum phase transition from a 
spin-liquid state to a long-range ordered antiferromagnet 
as td/t is reduced below a critical value, where the spin 
gap becomes zero^^ 



The key to the observation of quantum magnetism in 
our system is the presence of two different exchange en- 
ergy scales. Without dimerization, this situation also 
occurs for anisotropic simple cubic lattices with tun- 
nelling t along two axes and a stronger tunnelling tg 
along the third direction. In this case the symmetry be- 
tween neighbouring links is restored and the detected sin- 
glet and triplet fractions are the same for both merging 
configurations. We observe a clear population difference 
(Ps — Pto)/2 after loading a gas with entropies 5in below 
1.0 /cb into an anisotropic lattice, which increases to 4% 
for larger tunnelling ratios tg/t, see Fig. |4|l. The popula- 
tion difference is equal to the transverse spin correlator 
between neighbouring sites i and i + 1 along the strong 
tunnelling direction 



(5f5f+i)-(5f5f+i) = (ft-pt„)/2. 



(2) 



This quantity hence directly characterizes the fraction 
of atoms with antiferromagnetic ordering on neighbour- 
ing sites in the entire atomic cloud. Our observations 
also extend to weak lattices, where correction terms 
to the single-band Hubbard model become relevantPH. 
In this r egime a variety of magnetic phases have been 
predicte d^^ ' ^^ l 

The results can be compared to a second order high- 
temperature series expansion (Methods). We find good 
agreement in the regime of small anisotropics. For larger 
anisotropics the expansion breaks down, as the strong 
tunnelling and the temperature become comparable. In 
this regime we expect the temperature to lie between the 
large and small exchange scales J < k^T < J^. The sys- 
tem then behaves as an array of one-dimensional spin- 
ordered chains, without correlations between thenP^. 
Low-dimensional systems have been predicted to show 
enhanced nearest-neighbour correlations^^. 

For temperatures much larger than the strong ex- 
change energy the magnetic correlations should disap- 
pear. In Fig. (4)3 we study the dependence on the ini- 
tial entropy s^^ and find the correlations to vanish above 
2.5 /cb, where k^T ^ Jg is expected. 

Owing to the presence of the harmonic trap, most 
spin correlated atoms are located in the centre, where 
the filling is close to one particle per site. The density- 
normalized fraction of antiferromagnetic ordering is ob- 
tained when dividing by the fraction of atoms with two 
particles of arbitrary spin on adjacent sites. Under the as- 
sumption that all spin correlators {S^''^'^ S^J^{^) are equal 
- which applies if all symmetry breaking fields are much 
smaller than all other energy scales - the normalized spin 
correlator S can be directly obtained from the measure- 
ment of singlets and triplets (here n- is one for a single 
particle of any spin on site i and zero otherwise) 



S = 



Ps + 3pto 



(3) 



The normalized antiferromagnetic correlations along the 
strong tunnelling direction reach 25%, see inset Fig. |4^. 
This corresponds to approximately 5, 000 ordered atoms. 
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In this work, we have demonstrated the observation 
of quantum magnetism of repulsively interacting ul- 
tracold fermions in cubic lattices and investigated the 
dependence on temperature, lattice dimerization and 
anisotropy. The presented approach should also facil- 
itate the access to low temperature regimes for two- 
dimensional geometries. Furthermore, the tunable ge- 
ometry optical lattice allows the extension of our studies 
to spin-ladder systems, dimerized one-dimensional chains 
and zig-zag chains, where the interplay between quantum 
fluctuations and magnetic ordering plays a particularly 
important rol#^^. At even lower temperatures, the ex- 
istence of spin-liquids in honeycomb or triangular lattices 
could be investigated^. 



to the approximate expressions 

Ux CX 

ujy (X ^Jv^^ 0.81(FxVy/^x) + 0-78^z 
uj^ (X yjv^ + 0.81(1/xVy/1^x) + 1.24Fy. 

For %,xY,z = [3.7(1), 0.13(1), 9.8(3), 11.0(3)] ^R, as 
in Fig. l3^, the lattice contributes trapping frequen- 
cies of a;^**f727r = [62.7(9), 57(1), 54.3(3)] Hz. Addi- 
tionally, the optical dipole trap creates a trapping of 
^^TzV^^ = [30.7(2), 105.9(3), 34.6(2)] Hz, in ah data 
sets shown, except for Fig. pp andpk, where uj^^^'f^/27T = 
[28.1(2), 90.1(3), 31.6(2)] Hi:^ 



METHODS 



A. Preparation. 

After sympathetic cooling with ^^Rb, 2 x 10^ fermionic 
^^K atoms are transferred into an optical dipole trap op- 
erating at a wavelength of 826 nm. A balanced incoherent 
spin mixture of atoms in the Zeeman levels tuf = —9/2 
and —7/2 of the F = 9/2 hyperfine manifold is then pre- 
pared and evaporatively cooled^ ^. When taking data as 
a function of entropy, the gas is heated through inelastic 
losses by setting the magnetic bias field to a value close 
to the Feshbach resonance at 202.1 G. We measure the 
entropy per particle in the dipole trap sin using Fermi fits 
to the momentum distribution of the cloud after expan- 
sion. The field is finally increased to 221.4 G resulting 
in a scattering length of 106(1) a^. The optical lattice is 
subsequently turned on in 200 ms using a spline shaped 
laser-intensity ramp. 



B. Trapping potential. 

The red-detuned optical lattice is created by four retro- 
reflected laser beams at A = 1064 nm (refs [T4l and [T5]). 
This gives rise to a potential of the form 

V{x, y, z) = -V^cos^ikx + 9/2) - Fx co^^{kx) 
—Vy cos^ {ky) 

—2a^/VxVY cos{kx) cos{ky) cos (p 
-Vzcos\kz), (4) 

where V^^ Vx, Vy and Vz denote single-beam lattice 
depths (as calibrated using Raman-Nath diffraction on a 
^^Rb Bose-Einstein condensate), k = 27r/A and the mea- 
sured visibility of the interference pattern a is 0.90(5). 
The phase (p is stabilized to 0.00(3)7r, whilst is set to 
1.000(l)7r. Gravitation points along the y direction. 

The laser beams create an overall harmonic trapping 
potential which scales with the lattice depths according 



C. Detection lattice ramp. 

For measurements in the dimerized lattice, the lattice 
is ramped up in two steps. All beam intensities are lin- 
early increased over the course of 0.5 ms up to the point 
where = 30(1) £;r, = 40(1) £'r and ah other inten- 
sities in the xy plane are ramped such that the potential 
is not deformed. In a second linear ramp lasting 10 ms, 
the lattice is changed to a simple cubic geometry where 
%,x,y,z = [25(1), 0,30(1), 40(1)] i?R. 

Our observable locally projects onto the two-site eigen- 
states of individual dimers, which includes an admixture 
of double occupancies. We use an exact calculation of a 
two-site Hubbard model to estimate the adiabaticity of 
this ramp. The unnormalized singlet ground state of this 
system is given by 



4td 



(lt4)-i;,t)) + (in,o) + |o,t;)). 



There is hence a significant contribution of double oc- 
cupancy in the regime where U ~ td, as applies for the 
most strongly dimerized lattices investigated in this pa- 
per, whilst the contribution vanishes for the deep simple 
cubic lattice used for detection where U/t ^ 600. For 
the given ramp-times and including site-offsets due to 
the harmonic trapping potential, the probability of pop- 
ulating excited states during such a ramp remains below 
0.1% for all values of U and explored in the dimerized 
systems. 

In the anisotropic lattice, we directly ramp to l^c y z ~ 
[25(1), 30(1), 40(1)] £;r in 0.5 ms {Vx = Er through- 
out). Our detection method then corresponds to lo- 
cally projecting the wavefunction of the system onto 
(It^i) ~ \iA))/V2 on pairs of sites when measuring p^^ 
hence excluding any contributions from double occupan- 
cies. The probability of this projection is above 80% for 
all shown data and is higher for deeper lattices. For both 
lattice geometries, the triplet state stays unaffected by 
changes in U and t. 
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D. Singlet-triplet oscillations. 

For both systems, once the ramp to the deep sim- 
ple cubic lattice has been completed, double occupan- 
cies are removed via spin-changing collisions, which oc- 
cur after transferring atoms from the tuf = —7/2 to 
the rriF = —3/2 state^. We verify that this proce- 
dure removes all double occupancies but leaves singly 
occupied sites unaffected by measuring the number of 
double occupancies and the total number of atoms be- 
fore and after applying the removal sequence. A mag- 
netic field gradient causing spin- dependent energy offsets 
of A_9/2,_7/2//i = [1291(1), 1156(1)] Hz between neigh- 
bouring sites is then applied, giving rise to coherent os- 
cillations between singlets and triplets^^. Subsequently, 
pairs of adjacent lattice sites are adiabatically merged 
into one by linearly reducing to zero whilst increasing 
Vx to 25(1) £^R. The double occupancy in the merged 
lattice is then measured as in ref. JT, but taking into 
account an independently calibrated detection efficiency 
of 89(2)%. We verify that merged sites containing two 
atoms of opposite spin but in different bands are not 
detected as double occupancies by artificially creating a 
state containing large amounts of triplets but no singlets 
and measuring t^d- 

We apply a sinusoidal fit to the double occupancy 
where the frequency and phase are fixed and take into 
account the damping of the oscillations, which was cali- 
brated independently and is included by multiplying the 
amplitude by 1.16. A phase shift arises owing to a weak 
residual magnetic field gradient present during the dou- 
ble occupancy removal procedure, whereas the contribu- 
tion from switching the singlet-triplet oscillation gradient 
on and off is negligible. We confirm that the maximum 
of the oscillation corresponds to its starting point (and 
hence to the number of singlets) by merging the lattice 
immediately after it has been ramped to a deep simple 
cubic structure. 
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SUPPLEMENTARY INFORMATION 

A. Dimerized cubic lattice - series expansion 

We start with the definition of the homogeneous single- 
band Hubbard Hamiltonian in a dimerized cubic lattice 
with strong and weak links between nearest neighbours 
(i, j)_ and (i, j) respectively, see Fig. [sj 

Hd = Ho + He 

Ho = -td + h-c-) + ^ X] ^it% 

-/i^(ni^ + ni;) (6) 

Hc = -t Yl (4,A- + h.c.). 

We have split the Hamiltonian into the dimer part Ho and 
the coupling between dimers He. The on-site interaction 
energy is given by /7, the tunnelling matrix elements be- 
tween nearest neighbours by t and td and the chemical 
potential is parametrized with /i. The fermionic creation 
operator for an atom on the lattice site i is given by cj^, 
where a G {t, i} denotes the magnetic sublevel and h.c. 




Dimerized cubic lattice Anisotropic cubic lattice 



FIG. 5. Lattice geometries. Overview of the three- 
dimensional lattice geometries explored in the experiment, 
a, In the dimerized cubic lattice the tunnelling between neigh- 
bouring sites is increased to td on a chequerboard pattern in 
the xy plane as compared to the weaker tunnelling t. The 
pattern replicates along the z axis, b, In the anisotropic cu- 
bic lattice the tunnelling is increased to tg along the x axis, 
while remaining t in the other two directions. 
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is the Hermit ian conjugate. The particle number opera- 
Denoting the inverse 



tor is hi = hi^ + n^;, 
temperature with /3 = 1//cbT, the thermal average of an 
observable O then reads in the grand canonical potential 



(O) 



TT{Oe-P"'^} 
TT{e-l^"^} 



(7) 



We now treat the coupling Hamiltonian He = tT as a 
perturbation, which leads to an expansion of the above 
expression in powers of the dimensionless parameter I3t 
(ref. i29,in the main text). The expansion is expected to 
be close to the exact result in the regime t <C /cbT <C 
td, U. For the partition function up to second order (de- 
nominator in Eq. [t]) we find 



Z = Zn 



f^' Jo Jo 



dr2dr,{f\r,)f\r2)W (8) 



The expression for the numerator is analogous. The par- 
tition function of the unperturbed Hamiltonian is de- 
noted by Zo, whereas (...)o denotes the thermal average 
of the unperturbed Hamiltonian 



(f '(ri)f '(t2))o = Tr{exp(-/?iJ„)f '(Ti)f '(t2)}/^o 



T'(r) 



Te 



-tHo 



(9) 



B. Dimerized cubic lattice - observables 



As the expansion is up to second order in the tun- 
nel coupling, it is sufficient to evaluate all expressions in 
a two-dimer basis. Denoting the single dimer Hamilto- 
nian in the grand canonical potential with the un- 
perturbed partition function then reads 



(10) 



The evaluation of the second order terms is done in a 
double dimer basis |^]^,^^), where |^^^) and |^^) each 
denote one of the 16 possible eigenvectors of the first and 
second dimer link. This essentially leaves the evaluation 
of matrix elements of the following kind 



{^]M\f'{n)f'{T2)\^,^]) 



and 



(11) 



which can be computed either directly or numerically. 
For the singlet and triplet fraction the observable O takes 
the form of a projector for the 16 possible states on 
a dimer. The entropy and particle number per dimer 
are evaluated from the grand canonical potential l^d = 
—k^TlogZd of a single dimer. Fig. |6] shows the dimer 
singlet probability and the entropy per site versus filling 
calculated in second order for different temperatures and 
dimerizations. A comparison between the predictions of 
lowest order (atomic limit) and second order is shown in 

Fig. [3 
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FIG. 6. Dimerized lattice. High-temperature series pre- 
dictions up to second order for the homogeneous dimerized 
cubic lattice. The dependence on filling of the dimer sin- 
glet probability and of the entropy per site is shown. We set 
JJ jt — ^ and td/t = 10 or 2 and use different temperatures 
k^T/t. The entropy at half filling for large dimer izat ion is 
strongly reduced. 




10 12 14 



ksT/t 



FIG. 7. Higher order contributions. Comparison of the 
high-temperature series predictions in lowest and second order 
(dashed and solid line) for the dimer singlet probability and 
for the entropy per site in a homogeneous dimerized cubic 
lattice. The filling is set to one particle per site and the 
interaction to U/t = 5. The second order contributions are 
expected to be larger for lower temperatures and lead to a 
reduction of the singlet probability. 



C. Anisotropic cubic lattice 

Similar to the case of the dimerized lattice, we split 
the Hamiltonian for the homogeneous anisotropic cubic 
lattice into two parts 



Ha = Hu + Ht 



(12) 



Ht = -ts (4^C^- ^ + h.C.) - t (4aCj,a+h.C.). 



Notations are analogous to the previous section, see 
Fig. [5] The strong tunnelling between nearest neigh- 
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bours along the x direction is denoted with t^, whereas 
the weaker tunnehing along the other two axes is given 
by t. We treat the tunnelling Hamiltonian Ht as a per- 
turbation to the unperturbed part Hu^ which leads to 
an expansion of the partition function as in Eq. |8] in 
powers of ^ts and ^t. Density and entropy per site are 
then obtained from derivatives of the second order grand 
potential 



f3Q = -logZi 



w 



(13) 

Here Zi is the unperturbed single site partition function, 
(" = exp(/3/i) the fugacity and w = exp(— /3/7). 

The evaluation of the two correlators defined in the 
main text {SfSfj^i) and S is slightly more complicated, 
as it involves two neighbouring sites. However, the coef- 
ficients for these correlators have already been computed 
(ref. |30]in the main text). 



S 



1 



w 



1 



1 w;- 1 



(14) 



(15) 



D. Harmonic trap 

The effect of the harmonic trap is included in a lo- 
cal density approximation with a quadratically varying 
chemical potential 



m(0 = Mo - ^mcJ^(^)V^ 



(16) 



where uJ is the geometric mean of the trapping frequen- 
cies, /io the chemical potential in the centre of the trap 
and r the normalized distance of a site to the trap cen- 
tre. Any trap averaged observable 0*^^p is then obtained 
from integration of the contributions per site 0^°"^(/i) 



0*^^P= / W(9^°^(/i(r))dr. 



(17) 



Owing to the harmonic trap, the energy offset between 
neighbouring sites on the dimer links changes over the 
cloud size. The relative correction of this effect to all 
shown quantities was computed to be less than a few 
percent. 



E. Singlet and triplet fractions 

The fraction of atoms forming singlets and triplets Ps 
and pto obtained from an integration over the left 
merging sites, which are denoted by A (as described in 
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FIG. 8. Dimerized simple cubic lattice, a, Singlet and 
triplet fractions on the strong dimer links vs. initial entropy 
before loading into the lattice. The data shown here is used 
to compute the normalized imbalance X displayed in Fig. [3^. 
b, Singlet and triplet fractions as a function of dimerization 
td/t, corresponding to the measurements of Fig.jsJ). Vertical 
error bars denote the fit error from singlet-triplet oscillations 
consisting of 63 measurements, the errors in td/t stem from 
lattice calibration uncertainties and the errors in sin are the 
standard deviation of five measurements. 



the main text, see Fig. [2] 

ps = 2J2{p:)/n P!=\^i){n (18) 

Pto =2^(f>'«)/iV =|**o)(*^«|. (19) 

ieA 

Here and are the projection operators on the sin- 
glet and triplet states |^-) and |^*°) on neighbouring 
sites i and i + 1, (...) denotes the thermal average and 
N the total atom number. For the measurements in the 
anisotropic simple cubic lattice, the projection operators 
are related to the spin operators Si = 1/2 . 
where a = (cr^^, cr^, cr^) are the Pauli matrices 



4 



SjSi 



i+1 



2SfSf 



+1 



(20) 
(21) 



From this the equality of Eq. |2]in the main text is imme- 
diately obtained. In the dimerized lattice the projection 
operator on the triplet reads the same, whereas for the 
singlet the two-site system needs to be diagonalized. 

Fig. [8] shows the individual atomic fractions of singlets 
and triplets {ps and ptg) measured in the dimerized lat- 
tice, which are used to compute the normalized imbal- 
ance X. The data in Fig. [8^ corresponds to the entropy 
scan of Fig. [3^ in the main text, whilst the set of mea- 
surements for different tunnelling ratios in Fig. [SfD cor- 
responds to Fig. |3]3. We additionally measured the frac- 
tion of atoms on doubly occupied sites D in the lattice 
after freezing out the atomic motion but before apply- 
ing the cleaning procedure and inducing singlet-triplet 
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oscillations. For the entropy scan we find a nearly con- 
stant value of D = 0.026(5). For the measurements taken 
at different td/t the double occupancy ranges between 
0.066(4) and 0.29(2). 
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Fig. |9] shows Ps and pt^ for the measurements in the 
anisotropic cubic lattice, which are used to compute the 
spin correlators — {SfSf_j_i) — {S^ Si 



and 



S. Fig.|9^( 



responds to the measurements versus entropy displayed 
in Fig. [4]3, whilst the scan of the tunnelling anisotropy 
shown in Fig. [SJd corresponds to the data presented in 
Fig. [4k. For the entropy scan D lies between 0.11(6) and 
0.21(1), whilst for the scan of ts/t we measure double 
occupancies between 0.14(2) and 0.19(2). 



FIG. 9. Anisotropic simple cubic lattice, a, Singlet and 
triplet fractions as a function of initial entropy before loading 
into the lattice, corresponding to the measurements of Fig. 
b, Singlet and triplet fractions vs. tunnelling anisotropy ts/t, 
which are used to compute the spin correlators presented in 
Fig. [4^. Vertical error bars denote the fit error from singlet- 
triplet oscillations consisting of 63 measurements, the errors 
in ts/t stem from lattice calibration uncertainties and the er- 
rors in Sin are the standard deviation of five measurements. 



